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Within classical transport theory, an approach is presented to derive systematically 
effective collision terms and noise sources for effective kinetic equations of non- 
Abelian plasmas. This procedure amounts to the 'integrating-out' of fluctuations, 
and is applicable for in- and out-of-equilibrium situations. When applied to a 
hot non-Abelian plasma close to equilibrium we recover the collision integral and 
the noise source first obtained by Bodeker. It is also shown that our approach, 
and hence Bodeker's effective theory, is consistent with the fluctuation-dissipation 
theorem. 



1 Introduction 

Understanding the physics of a quark-gluon plasma at high temperature or 
density - which is expected to be detected in the up-coming heavy ion ex- 
periments at RHIC and LHC - requires reliable theoretical tools to describe 
transport phenomena in non-Abelian plasmas in- or out-of-equilibrium. Such 
methods are also of relevance for the physics of the early universe, if the baryon 
asymmetry could finally be understood within an electroweak framework. 

In this talk[] an approach is discussed which allows for 8,-systematic deriva- 
tion of effective kinetic equations for non-Abelian plasmasEru The interest in 
developing a^ formalism based on a classical transport theory for non-Abelian 
gauge fields □ is that main properties of a hot quantum plasma can already be 
understood within simple classical terms. Indeed, the soft non-Abelian gauge 
fields - those having a huge occupation number - can be treated as classical 
fields, while the hard gauge field modes can be treated as classical particles. 
This method, based on 'integrating-out' fluctuations about some mean fields, 
has been known since long for Abelian plasmas u, but a consistent extension 
to the non-Abelian case was longtime missing. 

We apply our formalism to a hot plasma close to equilibrium. While the 
HTL effective theory is recovered in the Vlasov approximation^, the inclusion 
of leading order corrections due to fluctuations results in Bodeker's effective 
theoryu For a more extensivejjliscussion of this and related effective theories we 
refer to the discussion in ref. Some further approaches to study quark-gluon 
dynamics are considered in ref. 
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2 Microscopic transport equations 

The starting point for a classical transport theory is to consider an ensemble 
of classical point particles carrying a non-Abelian charge Q a , where the colour 
index runs from a = 1 to TV 2 — 1 for a SU(iV) gauge group. These particles 
interact self-consistently amongst each others, that is, through the classical 
gauge fields created by the particles. Their classical equations of motion are 
known as the Wong equations 113, 

m^=r, m ^= 9 Q a FrP», m^ = - 9 r bc rAlQ c . (1) 

Here, denotes the gauge field, F* v = M A£ - d v A^ + gf abc A^A c u the cor- 
responding field strength and z = {x,p, Q){t) the world line of a particle. An 
ensemble of such particles can be described by their one-particle distribution 
function / ~ J2i J drS(z — Zi). Making use of (]l|), its kinetic equation isQ 

P» ~ 9f^KQ^ a - 9 Q a F^) f( X , P , Q) = . (2) 

This Boltzmann equation is collisionless, since df /dr = (Liouville's theorem). 
However, it contains effectively collisions due to the long range interactions 
as shall become clear in the sequel. Eq. @ is completed by the Yang-Mills 
equation D^F^" = J u , with the current J[f] obtained self-consistently from 
the particles. In this classical picture, / is a deterministic quantity once all 
initial conditions have been fixed. This is, however, not possible, and / has to 
be considered as a stochastic (fluctuating) quantity instead. 

3 Coarse graining 

Of particular interest for most physical applications are the macroscopic prop- 
erties of such a plasma of particles. In this case, a more appropriate quantity 
to study would be a coarse-grained distribution function /, whose dynamical 
equation is expected to be of the Boltzmann-Langevin type, 

v» (A* - gQaW) f = C[f] + c . (3) 

(We introduced the shorthands D^[A]f = [<9 M - gf abc QcA^ b d$]f , <9 M = d/dx^, 
df t = d/dp^ and = d/dQ a ). A coarse-grained distribution function obtains 
from the microscopic one after a suitable smearing-out over some characteristic 
volume, which, when applied to (^), results in a collision integral C[f] and 
a related source for stochastic noise (. Physically speaking, these terms are 
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expected because the coarse-graining integrates-out the modes within a coarse- 
graining volume. The interactions of such modes can result into additional 
effective interactions for the remaining long range modes. Also, the stochastic 
fluctuations of / are smeared-out, resulting in a noise term (. 

4 Effective transport equations 

An effective kinetic equation, like ((^), can be derived systematically within 
classical transport theory, starting with (|^). To that end, we perform an 
ensemble average in phase space which allows to split / and the gauge fields 
into their mean field part and a fluctuating part 

f(x,p,Q) = f(x,p, Q) + 5f(x,p, Q) , (4) 
A»{x) = A»(x) +a^(x) , (5) 

where / = (/) and A = (A). The mean value of the statistical fluctua- 
tions vanish by definition, (<5/) = and (a) = 0. This separation into mean 
fields and statistical, random fluctuations corresponds effectively to a split 
into long/short wavelength modes associated to the mean fields/fluctuations. 
Performing the ensemble average of (|J) and the Yang-Mills equations gives uU 

v» (A. - gQaP^d;) f = (v) + (0 • (6) 

D„F^ + (J%J = r . (7) 

In (||) and (Q) , we collected all terms quadratic or cubic in the fluctuations into 

v = gQ a V 1 " [(D»a v - D v a„) a + g^afa] %6f(x, p, Q) , (8) 
£ = g P »f abc Q c [d®al Sf(x, p, Q) + ga«ald»f(x, p, Q)] , (9) 
J£ ee g [f dbc D^a b ,^ c + f abc a b ,, {{D»a v - D v a») c + gf cde a^ e )] . (10) 

In addition, we need the dynamical equations for the fluctuations. They are 
derived from subtracting (g) from (||) as 

p" (D M - gQaF^d;) Sf = gQ a (D^ - D v a») a p^f 

+gp^ b ^f abc Q c dQf + V + £ - ( v + (If) 
(5V - D»{D y a»)) a + 2gf abc F^a c , u + J%£ - (J££) _ SJ a >» . (12) 

Eqs. (|TT| ) and (|fj) are the master equations for all higher order correlation 
functions of fluctuations. (This hierarchy of dynamical equations is similar 
to the BBGKY hierarchy.) All the equal-time correlators of Sf can be de= 
rived from the basic definition of the Gibbs ensemble average in phase spaceB 
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The most important correlator is the quadratic one (dfdf)t = tuf^h^ch can be 
expressed in terms of / and a two particle correlation functionHO 

Some few comments, (i) The set of dynamical eqs. - ( |l2| ) is equivalent to 
the original set of microscopic equations. It is exact in that no approximations 
have been performed so far, and can be applied to both in- or out-of-equilibrium 
situations, (n) The original equations as well as their split into mean gauge 
fields and fluctuations have been shown to be consistent with the non-Abelian 
mean field gauge symmetryaQ (Hi) In order to obtain effective equations for 
the mean fields only, one has to integrate-out the fluctuations. This amount 
to solving explicitly the dynamical equations of 5f and the gauge field fluctua- 
tions for given initial conditions and to express the correlator terms explicitly 
as functions of the mean quantities, making use of the basic equal-time cor- 
relators. In this light, our procedure can be seen as a derivation of collision 
integrals, (iv) The dynamical equations for the fluctuations are non-linear as 
they involve fluctuations up to cubic order. Solving them will require some 
approximations. For small gauge coupling an expansion in powers of g can be 
performed. If the fluctuations remain small within a coarse-graining volume, 
the second moment approximation can be used, which amount ultimately to 
a linearisation of the dynamics of the fluctuations^ These types of approxi- 
mations are systematic and consistent with the non-Abelian gauge symmetry 
of the mean fields, (v) In the Abelian limit (£) vanishes identically, while (77} 
reduces to the well-known Balescu-Lennard collision integralB (vi) In the limit 
where fluctuations are fully neglected, the set of equations are known as the 
(non-Abelian) Vlasov equationsQi3 



5 The hot non-Abelian plasma close to equilibrium 

We now turn to the specific example of a hot thermal plasma close to equi- 
librium. 'Hot' implies that particle masses can be neglected to leading order 
m(T) <C T, and that the gauge coupling, as a function of temperature, is very 
small g(T) <C 1. The non-Abelian plasma is characterised by two dimensionless 
parameters, the gauge coupling g and the plasma parameter e, which obtains 
as the ratio of the volume per particle over a Debye volume. The Debye (or 
screening) length describes the polarisation of the plasma. A small plasma 
parameter is mandatory for a kinetic description to be viable. Physically, this 
means that a large number of particles interact coherently within a Debye vol- 
ume (quasi-particle behaviour). For a classical plasma, e is an independent 
parameter related to the mean particle number and the gauge couplings For a 
quantum plasma, the mean particle number can not be fixed arbitrarily, and 
e scales as e ~ g 3 <C 1. This explains why a kinetic description, for quantum 
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plasmas, is intimately linked to the weak coupling limit. 

To simplify our analysis we shall perform some approximations, all of which 
can be understood as a systematic expansion in g. After ensemble averaging, 
the distribution function / is effectively coarse-grained over a Debye volume. 
Fluctuations of / within a Debye volume are parametrically suppressed by 
powers of g. We expand / about the equilibrium distribution function / Gq to 
leading order in g as 

f = n + gf (1) +5f. (13) 

Solving (JsJ) using ([l3|) for Sf — reproduces the HTL effective theoryi Two 
further approximations are employed beyond the HTL level, when Sf ^ 0. We 
shall discard cubic correlator terms in (^|), which are additionally suppressed 
by a power in g, in favour of quadratic ones. This corresponds to neglecting 
three-particle collisions in favour of binary ones. We also employ the second 
moment approximation - setting r\ — (rj) and £ = (£) in ( |TT| ) - which means 
that the effects of collisions are neglected for the dynamics of the fluctuations. 
All these approximations can be systematically improved to higher order. 
Let us consider colour excitations, described by the colour current density 

J£(x,v)=gvi> J dQ d Po d\p\p 2 6(p ) S(p 2 ) Q a f(x,p,Q) , (14) 

and = p^/po — (1, v), v = 1. The current J is obtained after an angle 
average over the directions of v as J{x) = J J(x, v). The colour mea- 
sure dQ is normalised J dQ = 1 and contains the Casimir constraints of the 
gauge group, such as J dQ£LQb — CiS a h with C2 the quadratic group Casimir 
(for more details, see ref. nO). The approximate dynamical equations for the 
fluctuations become 

SJ p ) a = -m 2 D v p v" (D„ao - A>a M ) a ~ ff/a6c« M a* J c ' p , (15) 
(D 2 a» - D»(Da)) a + 2gf abc Fg v a c , v = SJ% , (16) 

where the Debye mass rnn obtains from the equilibrium distribution function 
and the group representation of the particles. Solving for the fluctuations in 
the present approximations yields a closed expression for 5 J, and an iterative 
expansion in powers of the background fields for aaa The seeked-for dynamical 
equation for the mean quasi-particle colour excitations is 

v^D^j + mlv^F^o = C lin [J°] + C(z, v) , (17) 

where the linearised collision integral Cii n can now be evaluated explicitly, using 
the expressions for a and SJ. The Yang-Mills equation remains unchanged at 
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this order. For the collision integral, one finally obtains to linear order in the 
mean current, and to leading logarithmic order (LLO) in (lnl/g) -1 <C 1 

CW[#](»,v) = 9f abc (al(x)8J^(x,v))\ A=0t LLQ 

= -iJ d ^i^y)j° a { x ,v) (is) 

with the kernel J(v,v') = <5 2 (v - v') - |(v.v') 2 /Vl = (v.v') 2 . Notice that 
the collision integral is local in coordinate space, but non-local in the angle 
variables. The rate 7 = ^-hxl/g is (twice) the hard gluon damping rate. 
In obtaining ([l8|), we introduced an infra-red cut-off ~ grriD for the elsewise 
unscreened magnetic sector. This result has been first obtained by Bodeker in 
ref. 0. In addition, the source for stochastic noise C in (|l5|) can be identified as 

Ux,v) = 9fa bc al(x)SJ c >»(x,v)\ A=oJ=0 (19) 

(C a (x,v)C b (y,v')) =2 im 2 D TI(v,v')6 ab S(x~y) , (20) 

and the self-correlator (|2(]) has been evaluated from (|l^) to LLO. Notice that 
its strength is determined by the kernel of the dissipative process. Next, we 
show the close relationship to the fluctuation-dissipation theorem (FDT). 



6 Fluctuation-dissipation theorem 

It is well-known that dissipation in a (quasi-stationary) plasma is intimately 
linked to the fluctuations, a link which is given by the FDT. For the coarse- 
grained kinetic equation (^) one may ask how the spectral functions of £ and / 
have to be related to C[f] in order to satisfy the FDT. This question has been 
considered in ref.d. Within classical transport theory, the FDT is implemented 
in a straightforward waytil The pivotal element is the kinetic entropy S[f]. We 
consider small deviations of / from the equilibrium, / = / cq +A/. The entropy, 
stationary at equilibrium, is expanded to quadratic order in A/, S = S cq + AS. 
We define the thermodynamical force as F(z\ = —6(AS)/5(Af(z)), where 
z = (x, p, Q). Given C[f](z), it can be shownalill that 

is the noise self-correlator for ( compatible with the FDT. For the particular 
example studied above, (20) follows from inserting (|l^) into (|2l]). This proves 



that Bodeker's effective theory is compatible with the FDT. Furthermore, the 
correlator (A/A/)| t=t / can be derived along the same lines and agrees, to lead- 
ing order, with (8fSf)\t=t> from the Gibbs ensemble average. This guarantees 
that our formalism is consistent with FDT. 
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7 Discussion 



We have presented a systematic approach to derive effective transport equa- 
tions for non-Abclian plasmas. The procedure amounts to the 'intcgrating-out' 
of fluctuations about some mean fields. Most importantly, the formalism is ap- 
plicable for both in- and out-of-equilibrium situations, simply because the main 
statistical information as encoded in the equal-time correlators of fluctuations 
does not depend on the mean distribution function being thermal, or not. At 
the same time, the formalism - and certain approximations to it - respects 
the underlying non-Abelian mean field symmetry, which is at the basis for any 
reliable computation. 

For the close-to-equilibrium plasma, we have shown how Bodekers effec- 
tive theory emerges to leading logarithmic accuracy. The collision integral 
as well as the necessary noise source have been derived explicitely from the 
microscopic theory. In addition, it is established on general grounds that 
the formalism (and hence Bodeker's effective theory) is consistent with the 
fluctuation-dissipation theorem. 

As a final comment we point out that the effective theory is the same for a 
classical or a semi-classical/quantum plasma, differing only in the equilibrium 
distribution function (Maxwcll-Boltzmann vs. Bose-Einstein or Fermi-Dirac) , 
and hence in the corresponding value for the Debye mass. The sole 'quantum' 
effect which entered our computation resides in the non-classical statistics of 
the particles, which is all that is needed to correctly describe a hot quantum 
non-Abelian plasma close to equilibrium at the present order of accuracy. 
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